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Abstract 

Zamolodchikov's recursion relations are used to analyze the existence and approxi- 
mations to the classical conformal block in the case of four parabolic weights. Strong 
numerical evidence is found that the saddle point momenta arising in the classical limit 
of the DOZZ quantum Liouville theory are simply related to the geodesic length functions 
of the hyperbolic geometry on the 4-punctured Riemann sphere. Such relation provides 
new powerful methods for both numerical and analytical calculations of these functions. 
The consistency conditions for the factorization of the 4-point classical Liouville action 
in different channels arc numerically verified. The factorization yields efficient numerical 
methods to calculate the 4-point classical action and, by the Polyakov conjecture, the 
accessory parameters of the Fuchsian uniformization of the 4-punctured sphere. 
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1 Introduction 



A few years ago a considerable progress in the Liouville theory has been achieved [1]. The so- 
lution to the quantum theory based on the structure constants proposed by Otto and Dorn [2] 
and by A. and Al. Zamolodchikov [3] was completed by Ponsot and Teschner [4-6] and by 
Teschner [1,7]. Along with the techniques of calculating conformal blocks developed by Al. 
Zamolodchikov [8-10] the DOZZ theory provides explicit formulae for quantum correlators. 

On the other hand there exists so called geometric approach originally proposed by 
Polyakov [11] and further developed by Takhtajan [12-15] (see also [16-19]). In contrast 
to the operator formulation of the DOZZ theory the correlators of primary fields with elliptic 
and parabolic weights are expressed in terms path integral over conformal class of Rieman- 
nian metrics with prescribed singularities at the punctures. The underlying structure of this 
formulation is the classical hyperbolic geometry of the Riemann surface. 

Although the relation between these formulations is not yet completely understood [20- 
22] it is commonly believed that the quasiclassical limit of the DOZZ theory exists and is 
correctly described by the classical Liouville action of the geometric approach. This is for 
instance justified by explicit calculation of the classical limit of the DOZZ structure constants 
and the classical Liouville action for the Riemann sphere with three punctures [3, 23]. 

Some of the predictions derived from the path integral representation of the geometric 
approach can be rigorously proved and lead to deep geometrical results. This can be seen as 
an additional support for the correctness of the picture the geometric formulation provides 
for the semiclassical limit of the DOZZ theory. One of the results of this type is the so 
called Polyakov conjecture obtained as a classical limit of the Ward identity [24-29]. It 
states that the classical Liouville action is a generating function for the accessory parameters 
of the Fuchsian uniformization of the punctured sphere yielding an essentially new insight 
into this classical long standing problem. Its usefulness for solving the uniformization is 
however restricted by our ability to calculate the classical Liouville action for more than 
three singularities. 

The existence of the semiclassical limit of the Liouville correlation function with the 
projection on one intermediate conformal family implies a semiclassical limit of the BPZ 
quantum conformal block [30] with heavy weights A = Q 2 5, Aj = Q 2 5i, with S, Si = 0(1) in 
the following form: 



1+6Q 2 ,A 



A 3 A 2 
A 4 Ai 



(x) ~ exp{QV4§S](*)}- (1.1) 



The function fg 5 ^ ^ (x) is called the classical conformal block [3] or (a bit confusing) the 
"classical action" [9, 10]. 

The existence of the semiclassical limit (1.1) was first postulated in [9, 10] where it was 
pointed out that the classical block is related to a certain monodromy problem of a null vector 
decoupling equation in a similar way the classical Liouville action is related to the Fuchsian 
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uniformization. This relation was further used to derive the A — ► oo limit of the conformal 
block and its expansion in powers of the q variable. 

The 4-point function of the DOZZ theory can be defined as an integral of s-channel 
conformal blocks and the DOZZ couplings over the continuous spectrum of the theory. In 
the semiclassical limit the integrand can be expressed in terms of the 3-point classical Liouville 
action and the classical block, and the integral itself is dominated by the saddle point A s = 
Q 2 5 s {x). One thus gets the factorization [3] 



S {cl \6 4 ,6 3 ,S 2 ,S i; x) = 

S^(6 4 ,6 3 ,6 s (x)) + S^(6 s (x),6 2 ,6 1 ) - f 5s(x) 



(1.2) 



<53 ^2 
($4 (5i 



(x) ~ fs s (x) 



53 S 2 

54 Si 



where 8^(64,63,62, Si;x) is the classical action for the weights 61,62,63, 6 4 located at 0, x, 1 
and 00, respectively, and 8^(63,62,61) is the classical Liouville action for the weights <5i , <^2 , ^3 
at the locations 0, 1, 00. Since the semiclassical limit should be independent of the choice of 
the channel in the representation of the DOZZ 4-point function one gets the consistency 
conditions 



S (cl \6 4 , 63, 6 s (x)) + S^(6 s (z),6 2 , 6 1 ) - f Ss(x) [g g] ( x ) - f 6s(x) 
= (6 4 ,6 u 6 t (x)) + (5 t (x),6 2 ,6 3 ) 



S3 S2 

<5 4 81 



(x) 



(1.3) 



fs t (x) 



Si S2 

S4, 5 3 



(1 - x) - f St{x) 



Si 5 2 
5 4 5 3 



X 



= 26 2 log xx + (6 1 ,6 3 ,6 u (x)) + S {c V (6 U {x),6 2 ,6 4 ) 



fs u (x) 



5 3 62 
81 <5 4 



fs u (x) 



S3 S 2 
Si S4 



where the saddle weights 6t(x), 5 u (x) in the t- and n-channel are simply related to the s- 
channel saddle point: 

6 t (x) = 6 S (1 - x), 6 u (x) 



6 I - 



As the conditions (1.3) follow from the semiclassical limit of the bootstrap equations of the 
quantum DOZZ theory we shall call them the classical bootstrap equations. 

The question arises what is the geometric interpretation of the saddle point conformal 
weight 6 s (x). Let us recall that the classical solution describes a unique hyperbolic geometry 
with singularities at the locations of conformal weights. For elliptic, parabolic and hyperbolic 
weights one gets conical singularities, punctures and holes with geodesic boundaries, respec- 
tively [23, 29, 31]. In the latter case the (classical) conformal weight 6 is related to the length 
I of the corresponding geodesic by 

,i - 1 L | _ 1 (i.4) 



4 4 



2tt 



where the scale of the classical configuration is set by the condition R = — ^ imposed on the 
constant scalar curvature R. 



2 



Let us note that the relation (1.4) appears also in the context of the quantization of 
the Teichmiiller space ([32], [33] and references therein). With this relation assumed one 
can show that the Hilbert space arising in the quantization of the Teichmiiller space of a 
Riemann surface and the space of conformal blocks of the Liouville theory on this surface are 
isomorphic as representations of the mapping class group. 

In the case of 4 singularities at the standard locations 0, x, 1, oo there are three closed geo- 
desies T s , r^, T u separating the singular points into pairs (x, 0|1, oo), (x, 1|0, oo) and (x, oo|0, 1) 
respectively. Since the spectrum of DOZZ theory is hyperbolic the singularities corresponding 
to the saddle point weights Si(x) are geodesic holes. One may expect that these weights are 
related to the lengths £{ of the closed geodesies Tj in corresponding channels: 



If the formula above proved true it would provide a new powerful tool for calculating the 
geodesic length functions which play an important role in analyzing the structure of the 
moduli space of Riemann surfaces ( [34] , and references therein) . It would also pave a way for 
an explicit uniformization of at least 4-punctured sphere. 

Our aim in the present paper is a numerical verification of the three conjectures mentioned 
above: the asymptotic (1.1), the classical bootstrap equations (1.3), and the relation (1.5). 
Up to our knowledge no rigorous proof of any of these relations is known. The basic difficulty 
is the conformal block itself which except of some special cases is only known as a formal 
power series in the x variable. Since the coefficients are defined in terms of inverses of the 
Gram matrices of a Verma module a direct analysis is prohibitively difficult. 

A more efficient method based on a recurrence relation for the coefficients was developed 
by Al. B. Zamolodchikov [8]. The assumption crucial for this method is that the formal 
power series defining the conformal block converges. It is actually believed (and well justified 
by all special cases where the conformal block can be calculated explicitly) that the radius 
of convergence is 1. Another commonly accepted hypothesis, still waiting for its rigorous 
proof, states that the only singularities of the conformal block with respect to the z variable 
are branching points (in general of a transcendental kind) at 0, 1, and oo [36]. This implies 
that the conformal block is a single-valued analytic function on the universal covering of a 
3-punctured Riemann sphere and can be expressed by a power series convergent in the entire 
domain of its analyticity [10]. A recurrence relation for calculating coefficients of this so called 
g-expansion [10] provides an extremely efficient method for numerical analysis of conformal 
block and was applied for testing the conformal bootstrap equations [3, 37]. Let us note that 
both the x- and g-expansion recurrence methods involve formulae which successfully passed 
numerical tests but a formal proof of them is still lacking. 

The organization of the paper is as follows. In Sect. 2 we briefly present the classical 
Liouville action and its relation to the monodromy problem of Fuchsian uniformization. In 
Sect. 3 the classical conformal block is introduced and its relation to a certain monodromy 




(1.5) 
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problem of the null- vector decoupling equation is clarified. These sections contain a material 
which is basically known and were added mainly for completeness and to set the notation. 

In Sect. 4 the Zomolodchikov's recurrence method is applied to calculate the x- and q- 
expansion of the classical conformal block in the case of parabolic external weights up to 
7th and 16th power, respectively. It is also checked by symbolic calculation that up to these 
orders the asymptotic (1.1) is correct. 

In Sect. 5 the saddle point weights and momenta are defined and numerically calculated 
in same special cases in which the lengths of geodesies in each channel are exactly known. 
An excellent agreement is obtained using both the x- and (/-expansions of classical conformal 
block. It is also verified that the geodesic length function numerically calculated from the 
formula (1.5) satisfies the upper and lower bounds known to mathematicians [35]. The 
calculations of this section provide a strong evidence that the formula (1.5) is correct and 
along with the (/-expansion for classical conformal block give an extremely efficient method 
for numerical calculations of the geodesic length functions on the 4-punctured sphere. 

In Sect. 6 some numerical tests of the classical bootstrap are presented. It is shown that 
it is satisfied with high precision, improving with the number of terms taken into account 
in the (/-expansion of the conformal block. This allows for calculating the classical Liouville 
action for 4-puncture sphere and, by the Polyakov conjecture, the accessory parameters of 
the Fuchsian uniformization. Finally, Sect. 7 contains conclusions and discussion of possible 
extensions of the present work. 



2 Classical Liouville action and the Polyakov conjecture 

A conformal factor of the hyperbolic metric on the surface X parametrized with a complex 
coordinate z is a solution of the Liouville equation 

d z d- z <t>(z,z) = ^e^ (2.1) 

Given a genus g of X and a set of points z\, . . . ,z n removed from X this metric is determined 
by the singular behavior of 4> at Zj—s. In the present paper we shall consider the case of X 
being a punctured sphere (g = 0) and choose a complex coordinates on X in such a way that 
z n = oo. The existence and the uniqueness of the solution of the equation (2.1) on the sphere 
with the elliptic singularities, 

f -2(l-fj)k)g|2-Zji + 0(l) as z^ Zj , j = l,...,n-l, 
\ -2(l + £„)log|z| + 0(l) as z^oo, 

was proved by Picard [38,39] (see also [40] for a modern proof). The solution can be inter- 
preted as a conformal factor of the complete, hyperbolic metric on X = C \ {z\, . . . , z n -i} 
with the conical singularities of the opening angles < < 2ir at the punctures Zj. The 
solution is known to exist also in the case of parabolic singularities (corresponding to £j — > 0), 
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with the asymptotic behavior of the Liouville field of the form 

J -2\og\z-Zj\ -21og|log \z-Zj\\ +0(1) as z -> zj, 

= 1 01 II Ol II I II j ^ m ( 2 - 3 ) 

I — 2 log \z\ — 2 log | log |z|| + 0{1) as z — ► oo. 

The central objects in the geometric approach to the quantum Liouville theory are the 
partition functions on X : 

{X) = j Vfie-Q 2 ^, (2.4) 
M 

where M, is the space of conformal factors on X with the asymptotics (2.2) or (2.3), and the 
correlation functions of the energy-momentum tensor 

<T(«i) . . . T(u fc )T(«;i) . . . T(w{) X) = Jv<f> e ~ Q2s ^ T( Ul ) . . . T(u fc )T(«>i) . . . T(wi), (2.5) 



with 



-- (a u ^(u,u)) 2 + a^(n,u) 



(2.6) 



The singular nature of the Liouville field at the punctures requires regularizing terms in 
the Liouville action: 

S h [4>] = 

d 2 z [|cty| 2 + A*e*] + £ (1 - y |<fe| k 2 <A + (1 + £„) y \dz\ k-A (2.7) 

*e i=1 \z-zj\=e \z\=\ 

n-1 

-2tt J] (1 - O) 2 log e - 2tt (1 + £ n ) 2 loge, 
j'=i 

where X e = C \ j (UjU l z ~ z jl < e ) u (M > 7) } • The form of ( 2 - 7 ) is valid for parabolic 
singularities (with corresponding £j = 0) as well. 

One can check by perturbative calculations of the correlators (2.5) [12] that the central 
charge reads 

C =1 + 6Q 2 . (2.8) 

The transformation properties of (2.4) with respect to the global conformal transformations 
show [12] that the punctures behave like primary fields with the dimensions 

A? = A,- = ^ (1 - £ 2 ) . (2.9) 

As for fixed £j the dimensions scale like Q 2 , the punctures correspond to heavy fields of the 
operator approach [3]. 

In the classical limit Q 2 — > 00 with all 

* = ^ = ^ (2.10) 
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kept fixed we expect the path integral to be dominated by the classical action: 



(X) 



~ e 



-Q 2 S( cl >(<5 i; ^) 



(2.11) 



where the classical action S( cl \5i ; Zj) is the functional SlH (2.7) evaluated at the classical 
solution (p of (2.1) with the asymptotics (2.2) or (2.3). Similarly: 



(T(z)X)~T c \z) e ^ (cl) (^), 

where T cl (z) is the classical energy-momentum tensor. 
From (2.6) and (2.2) or (2.3) it follows that 

A, 



(2.12) 



T c \z) 
T c \z) 



and consequently 



A 

z- 

n-l 



T c \z)=Q 2 Y J 



,,\2 



for z — ► Zj, 
for z — ► oo, 



+ 



Combining now (2.11), (2.12) and (2.14) with the conformal Ward identity [30] 



(2.13) 



(2.14) 



n— 1 



(T(z)X) = J2 

3=1 



A, 



+ 



1 



d 



we get the relation 



{z — Zj) 2 z — Zj dzj 

dS^(5 l ;z l ) 
dzi 



(X), 



(2.15) 



(2.16) 



known as the Polyakov conjecture. 

It is amazing that this relation obtained by general heuristic path integral arguments 
turned out to provide an exact solution to a long standing problem of the so called accessory 
parameters of the Fuchsian uniformization of the punctured sphere. Indeed, it was rigorously 
proved [24-29] that the formula (2.16) yields the accessory parameters Cj for which the 
Fuchsian equation 

(2.17) 



<) 2 r(;)+± i T dL (zMz)=0 



admits a fundamental system of solutions with SU(1, 1) monodromies around all singularities. 
Note that if {xi{z),X2{z)} is such a system then the function <p(z,z) determined by the 
relation 



4 \wf 



/|2 



w(z) = 



Xi{z) 



(2.18) 



^(1-H 2 ) 2 ' y ~' X2{z)' 

satisfies (2.1) and (2.2) (or (2.3)). The SU(1, 1) monodromy condition is then equivalent to 
the existence of the well defined hyperbolic metric on X. 
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From (2.13) it follows that c,- satisfy the relations 



■j 

n—l n—l 



E c , =0 > J2(6 j + z j c j ) = 6 n . (2.19) 

3=1 3=1 

In the case of three singularities the only two accessory parameters are completely determined 
by (2.19). Hence one can solve the equation (2.17), find the classical Liouville field and 
calculate the classical action. The result is [3]: 

3^(63,62,5,) = S^(5 3 ,5 2 ,5 1 ;^,l,0) = (2.20) 

\(l - 6 - 6 - 6) log^+ F (^IT 1 + *2§ + *sf ) - £ + const, 



where 



X 

F(x) = J 



dy log- 



r(y) 



r(i-y)' 

1/2 

For the elliptic and parabolic singularities one has (5 ^ \. As we shall see below, it is useful 
to know the classical Liouville action also for 5 > \, what corresponds to the hyperbolic 
singularities. The relevant construction of the metric and the classical action was given in 
[23] . If we write 

6 3 = \+p 2 , p £ M, (2.21) 
then in the case of two elliptic/parabolic and one hyperbolic singularity 

sw (83,82,81) = ^(i-ei-6)iog^ + E i? (^ i + ^| + ^3p) 

<72,<T3 = ± 

2 

- J2 F fe) + H ( 2i P) + Ap\ + const, (2.22) 
3=1 

with 

dy log 



H( X ) = Jdy log^. 



3 Classical conformal block 

The partition function (2.4) corresponds in the operator formulation to the correlation func- 
tion of the primary fields V aj (zj, Zj), 



(X) = (V an (^, 00) . . . V ai ( Zl ,zi)), (3-1) 

where 

Aj = aj (Q - ay), ay = ® (1 + Zj) ■ 
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The DOZZ 4-point correlation function with the standard locations z = 0, x, 1, oo is expressed 
as an integral over the continuous spectrum 



V aA (oo, oo)y« 3 (1, l)V a2 (x, x)V ai (0, 0) 



(3.2) 



da C(a 4 , a%, a)C(Q — a, a 2 , oi\) 



1+6Q 2 ,A 



A 3 A 2 
A 4 Ai 



(s) 



Let 



1a,a = £(|£a,i> ® |€a,i»«€a,i| ® <£a,: 



be an operator that projects onto the space spanned by the states form the conformal family 
with the highest weight A. The correlation function with the 1a, A insertion factorizes into 
the product of the holomorphic and anti-holomorphic factors, 



y 4 (oo,oo)V3(l,l)lA,A^ 2 (x,x)Vi(0,0)^ = 
C(a4,a 3 ,a)C(Q - a, a 2 , «i) ^i+6Q 2 ,A 



(3.3) 



A 3 A 2 
A 4 Ai 



1+6Q 2 ,A 



A 3 A 2 
A 4 Ai 



(x). 



Assuming a path integral representation of the l.h.s. one should expect in the limit Q — > oo, 
with all the weights being heavy A, A, ~ Q 2 , the following asymptotic behavior 



F 4 (oo,oo)V3(l,l)lA,A^(x,x)y 1 (0,0)\ ~ e -Q 2 s^H^). 



(3.4) 



On the other hand one can calculate this limit for the DOZZ coupling constants [3, 23] 
obtaining 

C(a A ,a 3 ,a)C(Q - a,a 2 , ai ) ~ e -Q 2 (s«(S 4 ,5 3 ,5)+S«(5,5 2 A)) . (3 . 5) 



It follows that the conformal block should have the Q — ► oo asymptotic (1.1) so that 



S^(5 t ,x;5) = S( cl \8 4 ,S 3 ,S) + S^(S,6 2 ,S 1 ) - [g (x) - f s 



$3 <>2 

84 Si 



(x). 



(3.6) 



It should be stressed that the asymptotic behavior (1.1) is a nontrivial statement on the 
(quantum) conformal block. Although there is no proof of this property yet it seems to be 
well justified by sample numerical calculations, as well as by its consequences. We shall briefly 
describe two of them. 

The first one is the classical bootstrap mentioned in the introduction. In the semiclassical 
limit the l.h.s of the formula (3.2) takes the form e~® 2s °' O^^fo^i;^ w h ere 

S (cl) (<5 4 , S 3 ,S 2 ,S 1 ;x) = S^(5 4 , S 3 , S 2 , J i; 00, 1, x, 0) . 

The r.h.s. of (3.2) is in this limit determined by the saddle point approximation 

00 
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where 5 S = i + p 2 and the saddle point Liouville momentum p s is determined by 



dp 



S( cl \6 u x;l+p% =pa = 0. 



(3.7) 



One thus gets the relation (1.2) first obtained in [3] and the classical bootstrap (1.3) as its 
consistency condition. 

The second implication of the asymptotic (1.1) is the relation of the classical block to cer- 
tain monodromy problem, which in fact proved to be essential for developing a very effective 
recursive method for calculating the conformal block itself (and therefore also its classical 
asymptotic) [8, 10,36]. Consider the null vectors on the second level of the Verma module, 
given by 



|x±> = [L-2 ~ 



2(2A ± + 1 



L 



|A±), 



(3.8) 



where |A±) are the highest weight states with (for c > 25) 

A± = ^(5-c±V(c-l)(c-25)). 
This expression simplifies if we use convenient parametrization of the central charge 



c=l + 6Q 2 , Q = l + b, 



(3.9) 



so that 



A_ 



- - V, 

2 4 ' 



A. 



1 3 

2 ~ 462' 



It follows from (3.8) that the correlators 



(x±(z)X) A 



def 



y 4 (oo, 00)1/3(1, 1)1a,aX± {z)V 2 (x, s)Vi(0, 0) 



(3.10) 



where x±(z) are the null fields corresponding to |x±), satisfy the null-vector decoupling 
equation: 

r 2 



dz 2 



+ 7± 



a 



1 — z I dz 



7± 



Ai 



+ 



A, 



+ 



A, 



( X ±(z)X) A = 
A± 



(3.11) 



+ 



+ 



x(l — x) 







with 



and 



(z — x) 2 (1 — z) 2 z(l — z) z(z — x)(l — z) dx 
A± = Ai + A 2 + A 3 + A± - A 4 



<X±(^>A. 



7± = -(2A± + l) 



7+ = ~b 2 , 7- = -p- 



For Q —>■ 00 we have either b —>■ or b —>■ 00. To fix the notation we shall concentrate on the 
first possibility. For 6^0 the operator with the weight A + remains "light" (A + = 0(1)) 
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and its presence in the correlation function has no influence on the classical solution of the 
field equations 4 . Consequently, for b — > 



(x + (z)X) A ~ x c \z) e '^{^H^s 3 ,s )+ s^(s,s 2 A)~fs[Z Z\ i?] <*>) 



and we get from (3.11): 
where 



<5i <5 2 



z 2 (z — x) 2 (1 — z) 



^3 | ^1+^2 + ^3-^4 | X(1-X)C(X) 



z(l — z) z(z — x)(l — z) 



d 



(3.12) 



X c \z) = 0, 

(3.13) 

(3.14) 



The accessory parameter C{x) can be determined from the following monodromy problem. 
Consider the correlation function 



G(z) 



dcf 



F 4 (oo, 00)^3(1, l)x+(*)Vk(0,0) 



(3.15) 



where Va is the primary field corresponding to the highest weight state £a- The null vector 
decoupling equation for this correlator reads 



b 2 dz 2 + 



1\ d A 
- - ^ + ^7 + 



1 — z z J dz z 2 (1 — z) 2 z(l — z) 

Substituting into this equation the most singular term in the OPE 

X+(z)Va(0,0) ~z k Va,(0,0) 

we get 



A3 + A + V A r A4 l^)=0. (3.16) 



k(k-1) + 6 2 (A-k) = 0. 



(3.17) 



In the limit 6^0: 



and (3.17) transforms to 



A = ^ ( i-a 



46 2 



(i-e 2 ) + o(i) 



For z — > e 2 ™ z : 













--I 











1-g 

Z 2 

i±i 

Z 2 



and (minus) the trace of the monodromy matrix (an invariant with respect to the choice of 
the basis in the space of solutions of (3.16)) is equal to 

2 cos 7r£. 



'This is reflected by the fact that limit b — > of the equation (3.11) exists only in the "+" case. 
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Note that we get the same monodromy invariant if we replace in (3.15) the primary 
field Va with any of its descendants. Note also that (in the classical limit) the monodromy 
invariant of the two independent solutions of (3.16) for a curve encircling is by construction 
equal to the monodromy invariant for a basis in the space of solutions of (3.13) along a curve 
encircling both and x (all point on this curve can be taken "to the left" of the operator 
1a)- This condition fixes the accessory parameter C{x) that appears on the r.h.s. of (3.13). 



4 Zamolodchikov's recursion methods 



The BPZ 4-point conformal block is defined [30] as a formal power series 5 



cA 



A 3 A 2 
A 4 Ai 



(*) = * A - A -Mi + J>^ c ; 



n=l 



A 3 A 2 
A 4 Ai 



(4.1) 



Studying the analytic structure of its coefficients AI. Zamolodchikov derived the recursion 
relation [8]: 



-r-n 

^cA 



A 3 A 2 
A 4 Ai 



9A 



A 3 A 2 
A 4 Ai 



prs 



A 3 A 2 
A 4 Ai 



C-C rs (A) c rs (A),A+r S 



A 3 A 2 
A 4 Ai 



(4.2) 



r>2 s>l 
n > rs >2 



where 



c rs (A) = 13-6 [T rs (A) + 



1 



and g£ 



A 3 A 2 
A 4 A 1 



T rs (A)J' 

T /^n _ rs-l+2A+V(r-s) 2 +4 (rs-1) A+4A 2 

are coefficients of the expansion of the hypergeometric function, 

oo 

2 -Pi (A + A 2 — Ai, A + A 3 — A 4 , 2A, x) = gR A J 

n=0 



An exact form of the coefficients R£ 



A 3 A 2 
A 4 Ai 



(see Appendix) was partially derived and par- 
tially guessed in [8]. Although no proof of this form exists it is well justified by numerical 
calculations. 

Once the expansion (4.1) is known one can calculate the coefficients of the power expansion 
of the classical conformal block 



fs 



<5 3 &2 
84 Si 



{x) = ( 6 - 6l - §2) i ogx+ j2x n f s n 



n=l 



S3 82 
84 Si 



(4.3) 



5 It is believed that the series (4.1) converges on the unit disc but up to our knowledge there is no rigorous 
proof of this fact yet. 
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directly from the asymptotic (1.1) 



S3 $2 

Si Si 



n=l 



A 3 A 2 
A 4 Ai 



(4.4) 



n=l 



where on the r.h.s. one first expand the logarithm into a power series and then the limit is 
taken for each term separately. 

In the present work we were interested in a special case of all parabolic external weights 
Aj = the central charge c = 1 + 6Q 2 , and the intermediate weight A = Q 2 Q + p 2 ) 
parameterized by the Liouville momenta p. We have checked by symbolic calculations that 
up to m = 7 the limits in (4.4) exist yielding explicit formulae for the first seven coefficients 
of the expansion (4.3). Up to the first five terms: 



f(p,x) = 

+ 
+ 
+ 

+ 



i i 

4 4 

1 _1 

4 4 



(x) 

(p 2 -i)logx+(i + 4 

f_9_, 13 p 2 , 1 

^ 128 "•" 64 "•" 1024 (1+p 2 

fl9_, 23 p 2 , 1 
^384 "•" 192 "•" 1024 (1+p 2 ) 

( 1257 , 2701 p 2 
\ 32768 



(4.5) 



32768 



8388608 (1+p 2 ) 2 + 8388608 (1+p 2 ) "•" 8388608 (4+p 2 ) 

1 



2097152 (1+p 2 ) 
7439 



+ 



81 



) 



( 2573 , 5057 p 2 
V 81920 " r 81920 



+ 



1048576 (1+p 2 ) 3 

j 3343 j 81 \ x 5 , 

2f 4194304 (1+p 2 ) 4194304 (4+p 2 ) ^ X " r 



4194304 (1+p 2 

The limitation of the formulae (4.1) and (4.3) is that the power series involved are supposed 
to converge only for \x\ < 1. 

A more convenient representation of the conformal block was developed by Al. Zamolod- 
chikov [10] who proposed to regard it as a function of the variable 



q(x) 



K(l-x) 



dt 



V(l-t 2 )(l-xti) 



The map C \ {0, 1, cc} 3 x — > q(x) € D yields a uniformization of the 3-punctured sphere 
by the Poincare disc D. If the points 0, l,oo are the only singular points of the conformal 
block, then the block is a single valued function on D and the series in q variable converges 
uniformly on each subset {q : \q\ < e~ e < 1}. It was shown in [10] that the conformal block 
can be expressed as 



cA 



A 3 A 2 
A 4 Ai 



(x) = x~ 



(1-X) 



^-Ai-A 3 



x e 3 ( g y 



-4(Ai+A 2 +A a +A 4 ) 
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= x^" Al ~ A2 (l -x) 



A 3 A 2 
A 4 Ai 



U) 



-A1-A3 



^-2(Ai+A 2 +A 3 +A 4 ) 



X (16qr- — H cA 



A 3 A 2 
A 4 Ai 



where 



'c,A 



A 3 A 2 
A 4 Ai 



( q ) = l + J2(^q) n H^ A 



A 3 A 2 
A 4 Ai 



n=l 



The coefficients in (4.7) are uniquely determined by the recursion relation: 



Tjn 



A3 A 2 
A 4 Ai 



E AT 



A3 A 2 
A 4 Ai 



H 1 



A-A rs (c) c,A rs (c)+rs 



A3 A 2 
A 4 Ai 



(4.6) 



(4.7) 



n > 0, 



(4.8) 



r>l s>l 
n > rs >1 



where for the central charge parameterized as c = 1 + 6 (b + : 

A rs (c) = = — , - -If + - — + 



1-r 2 ,, 1 — rs 1 — s 2 1 
— — 6 2 + — 



and i?r s 



A 3 A 2 
A 4 Ai 



are related to the coefficients R A 



A 3 A 2 
A 4 Ai 



4 6 2 ' 



and their explicit form is known 



[8] (see Appendix). There are some important advantages of the formulae above. First of 
all one gets a series which is supposed to converge on the whole domain of analyticity of the 
conformal block. Secondly if we choose the parametrization (3.9) of the central charge the 
formulae does not contain square roots which simplifies symbolic calculations a lot. 

The representation (4.6) and the asymptotic (1.1) imply the following representation for 
the classical conformal block: 



fs 



S3S2 
84 <5i 



(*) = (l-S 1 -S 2 )logx + {\-S 1 -S 3 )\og{l-x) 
+ (l-2(5 1 + S 2 + S 3 + S 4 ))log(lK(x)) 



(4.9) 



+ (.5 -\) log 16 -(*-i) 7 rffi 7 -^ + /,; 



K{x) 



&3 <5 2 

84 81 



(?)■ 



where 



hs 



83 8 2 

84 81 



(q) = J2(^q) n K 



n=l 



83 82 

84 81 



(4.10) 



The coefficients in (4.10) can be calculated step by step from the asymptotic (1.1) 



E( 16 ^) nh s 



$3 <>2 
84 81 



1 



hm bg l + ^(16g)^ A 



n=l 



Q 2 ^oo Q 



A 3 A 2 
A 4 Ai 



(4.11) 



n=l 



Using this formula in the special case of all parabolic external weights we have checked the 
asymptotic (1.1) up to the terms q 16 . We have also calculated the coefficients of (4.10) up to 
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this order. The first few terms of the series (4.10) read: 



h(p,q) = hi 2 



+ 
+ 

+ 



i i 

4 4 

1 1 

4 4 



(q) 



1 + 



1 „2 

4(l+p 2 ) 9 



-1 , 

32(l+p 2 ) ;i ^ 



3 , 15 , 81 \ n 4 

128 (1+p 2 ) 2 ^ 128 (1+p 2 ) ^ 128 (4+p 2 ) ^ « 



(g6(l+; 



59 



(1+p 2 ) 5 384 (1+p 2 ) 4 256 (1+p 2 ) 3 + 2048 (1+p 2 ) 2 

1 661 _ 9 

16384 (1+p 2 ) 128 (4+p 2 ) 

-5 , 35 



<>(i.l 9 1 16875 ^ 6 

^ 16384 (9+p 2 ) / y 



1024 (1+p 2 ) 
1725 



+ 



4096 (1+p 2 )' 
1653 



+ 



141 

8192 (l+p 2 ) J 
_|_ 6279 



65536 (1+p 2 ) 4 65536 (1+p 2 ) 3 262144 (1+p 2 ) 2 

I 54913 6561 , 19683 

"^2097152 (1+p 2 ) 16384 (4+p 2 ) 3 262144 (4+p 2 ) 2 

, 332239 50625 , 1500625 \ i 

^ 1048576 (4+p 2 ) 2097152 (9+p 2 ) 1048576 (16+p 2 ) ) q 



+ ... 



(4.12) 



5 Geodesic length function 

According to the uniformization theorem the punctured sphere X is conformally equivalent 
to M./G, where EI is the upper half plane endowed with the Poincare hyperbolic metric and 
G is the Fuchsian group uniquely (up to conjugation in PSL(2,M)) determined by X and 
isomorphic to its fundamental group. 

The group G is generated by Tj G PSL(2, R), i = 1, . . . , n. It is possible to chose them in 
such a way that 

TiT 2 ...T n =I. 

If all the punctures correspond to parabolic singularities, then 

|TrTi| = 2. 

For each pair of punctures, say Z{ and Zj, there exists a unique closed geodesies of the 
hyperbolic metrics on X, separating Zi and Zj from the remaining singularities. Its length 
i (■fij) can be determined from the relation 

2 cosh = \Tr TiTj | . (5.1) 

Setting for the four punctured sphere (by an appropriate global conformal transformation) 
Z\ = 0, Z2 = x, z% = 1 and z 4 = co we have in the notation from the previous sections 

£ (712) = £ (7s) , £ (723) = £ (it) , £ (724) = £ (7«) • 
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In the case of locations x = g , \ , e *3 the group G is explicitly known [35] 



e 3 



Ti 



T 2 



T 3 



3 -8 
2 -5 



T 4 



Taking into account the locations obtained from x = |, 3 by SL(2, C) transformations 
preserving the set {0, 1, oo} one gets 



(5.2) 





1 18 9 Q Q 

9 8 9 8 ° y 


5 -1 2 


7ri i iri 

e 3 e 3 


cosh ^ 


2 2 5 8 5 8 


3 3 7 


7 7 
2 2 



All the cases listed above concern the metrics on M./G, induced by the hyperbolic Poincare 
metric on M with the scalar curvature —2. This corresponds to \i = 4 in the relation (1.5) 
which in terms of the s-channel saddle point Liouville momentum (3.7) reads 



Kls(x)) = 4irp s (x) 



(5.3) 



Taking into account the explicit formula (2.22) for the 3-point classical action one can write 
the saddle point equation (3.7) in the form: 

& T(l + 2ip)T 2 (±-ip) dp^ ' v ' 

The precision of a numerical solution to this equation is in practice determined only by the 
precision of an approximation to the classical conformal block f(p,x) at hand. 



For the locations inside the unit disc x 



l l 



g , 2 > § one can consider approximation 



given by the first n terms of the x-expansion (4.5). We have done the numerical calculations 
for n = 0, 1, . . . , 7. The deviations from the exact values measured by the differences 

A(x) = cosh [i.C(7 s (x))] — cosh [2irp s {x)} 

are presented in Tab.l. They provide a very good confirmation of the formula (5.3). 



A much better precision can be achieved if we use the (/-expansion of f(p,x) (4.12). 
The results of numerical calculations with the approximations of f(x,p) up to the terms 
q 2k , k = 0, 1, . . . , 8 are presented in Tab. 2. The agreement with the conjectured exact formula 
(5.3) is perfect. 
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n 




A(-|) 


A(|) 


A(f) 





1.8 x 1(T 2 


- 1.9 x 1(T 2 


2.8 x KT 1 


1.72 


1 


8.5 x 1CT 4 


9.4 x 1CT 4 


7.0 x 10~ 2 


1.05 


2 


5.6 x 1(T 5 


- 6.8 x 1(T 5 


2.2 x 1(T 2 


0.72 


3 


4.3 x 1(T 6 


5.8 x 10~ 6 


8.1 x 1(T 3 


0.53 


4 


3.6 x 10~ 7 


- 5.4 x 1(T 7 


3.1 x 10~ 3 


0.40 


5 


3.2 x 1(T 8 


5.3 x 1(T 8 


1.2 x 1(T 3 


0.31 


6 


2.9 x 1(T 9 


- 5.4 x 10~ 9 


5.3 x 10~ 5 


0.24 


7 


2.7 x 10~ 10 


5.7 x 1(T 10 


2.3 x 1(T 5 


0.19 



Tab.l 



2k 


A(|) 


AH) 


A(§) 


A(|) 





- 4.0 x 10~ 6 


- 4.0 x 10~ 6 


- 3.9 x 10~ 4 


- 1.1 x 10~ 2 


2 


- 1.0 x 10~ 1() 


- 1.0 x 10~ 10 


- 3.7 x 10~ 7 


- 1.0 x 10~ 4 


4 


- 1.8 x 10~ 15 


1.8 x 10~ 15 


- 5.4 x 10 _n 


- 2.3 x 10~ 7 


6 


- 1.3 x 10~ 15 


1.8 x 10" 15 


- 2.5 x 10 -13 


- 7.2 x 10 -9 


8 


- 1.3 x 10~ 15 


1.8 x 10~ 15 


1.8 x 10~ 15 


9.9 x lO" 11 


10 


- 1.3 x 10~ 15 


1.8 x 10~ 15 


1.3 x 10~ 15 


- 2.0 x 10~ 12 


12 


- 1.3 x 10~ 15 


1.8 x 10~ 15 


1.3 x 10" 15 


3.0 x 10~ 14 


14 


- 1.3 x 10~ 15 


1.8 x 10~ 15 


1.3 x 10~ 15 


4.4 x 10" 15 


16 


- 1.3 x 10~ 15 


1.8 x 10~ 15 


1.3 x 10~ 15 


1.5 x 10~ 14 



Tab. 2 



The advantage of the g-expansion is that it converges rapidly for all locations x on the 
complex plane except tiny neighborhoods of x = 0, l,oo. One can thus verify the formula 
(5.3) for all the cases listed in (5.2). Using the (/-expansion up to the terms q 16 we get the 
deviations from the exact values (5.2) closed to the precision of our numerical calculations 
set to be slightly less than 16 decimal digits (Tab. 3). 

As another verification of the formula (5.3) one can compare the numerical results with 
the analytic bounds for the geodesic length function for the locations x such that the geodesic 
7 S (x) is contained in the unit disc [35]: 

2 ~ < t>(x) < 7 . ; " (5-5) 



log (256IXI- 1 + 136) sv ' log ^ i6^T3^ j 



isyii-s| 

M 



41og^ 

We present sample numerical calculations of the geodesic length function along the rays 

x(r) = re i * k , fc = 0,l,2,3,4. (5.6) 



16 



X 


A(x) 


X 


A(x) 


1 

9 


- 1.3 x 1(T 15 


1 

2 


1.3 x 1(T 15 


1 

8 


i.8 x i(n 15 


-1 


4.4 x 10" 15 


8 
9 


i.5 x i(n 14 


2 


- 7.1 x 10 _1B 


9 
8 


1.4 x 1CT 13 


e _ "3" 


2.7 x 10 _1B 


-8 


- 8.9 x 1(T 15 




2.7 x 10" 1B 


9 


i.4 x i(n 13 



Tab.3 



The results for the (/-expansion up to the term q w are plotted on Fig.l for two different ranges 
of r. 

Let us close this section with the remark that it is straightforward to work out an analytic 
approximation for the saddle point momentum p s and, consequently, for the geodesic length 
£(j s ) f° r small x. Indeed, with the help of the Lagrange duplication formula 

T(2z) = (27r)-h 2z -h(z)T(z + ^) 

eq. (5.4) can be rewritten as 

-x + i6p.o g 2 + » log = (5.7) 

Performing a series expansion of the l.h.s. according to the formula 

r < 1+ «) I h * /' 

where 7 is the Euler-Marcheroni constant, we can rewrite (5.7) in the form 

-TT + 16plog2 + 8p 3 f; (-l) fc - 1 (l 2 -4 fc )C(2fe + 1) ^-2 = R |_ /(pj x) _ (5g) 
fc=i 

Using the explicit form of /(p, x) one can express a solution of (5.8) in the form of an expansion 
in 5ft x and Xo ^ xS For instance, keeping in (5.8) terms up to p 3 and x 3 we get 

TT 



— log 


; xx + 16 log 2 


-5ftx 




_205_Sft T 3 
1536 JLJ/ 




8C(3) + 


256 


x 2 + 5ft x 3 ) 




— log 


; xx + 16 log 2 


-5ftx 

7T 


_ 207«? T 2 _ 
512 JL x 


_205_<» r 3 
1536 JLJ/ 



log xx + 16 log 2 - 5ft x - fgJftx 2 - T^ 3 ^ 3 



p 3 + 
+ 



6 tlx 4 



log XX / ' log XX 



8vr 3 C(3) + |J(5ftx 2 + 5ftx 3 ) /7 1 V 5ftx 4 



(-logxx + 161og2-5ftx- gg5ftx 2 - ^5ftx 3 ) ^Vlogxxy logxx 
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2x10" 7 4x10" 7 6x10" 7 8x10" 7 1x10" 6 

Fig.l: Numerical calculations of the s-channel geodesic length function of the hyperbolic metric with the 
constant scalar curvature R — —2 along the rays (5.6) and the analytic bounds (5.5). 

6 Classical bootstrap 

The (/-expansion yields an extremely efficient method of numerical calculation of the classical 
conformal block and the saddle point momenta for all locations i / 0, l,oo. The precision 
of the s-channel calculations certainly worsens near x = l,x = oo singularities. Still the 
range of rapid convergence of the g-expansion is huge enough to provide a reasonable testing 
ground for the classical bootstrap equations (1.3). For instance in the region \x\ < 4.6 x 10 5 
with the small disc |1 — x\ < 0.0003 around x = 1 removed one has \q(x)\ < ^. 

In the present paper we are interested in checking the classical bootstrap equations in 
the simplest case of four parabolic singularities. The s-channel factorization of the 4-point 
classical action (1.2) yields the expression 

S,(x) = <?( cl )(!,i,i,i;x) 

= 2S 3 ( Ps (x))-2Uf(p s (x),x) 
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where f(p,x) is the classical conformal block (4.5), p s {x) is the saddle point momentum in 
the s-channel given by (5.4), and £3 (p) is the 3-point classical action (2.22) for two parabolic 
and one hyperbolic weight 6 : 

= 4F (i + ip) + H(2ip) + np, 

The efficiency of numerical methods provides a considerable freedom in testing the clas- 
sical bootstrap equations. As an example we have chosen the calculation of the relative 
deviation from these equations measured by the functions: 



A su (x) 



Sa(1 



Sa 



S^x) 

— log \x\ 



S 4 (x) 



S 4 (x) 



The largest deviation should be expected around the locations x = 0,1, 00. The behavior 
of the functions A st (x), A su (x) in these regions is well represented by their values along the 
real axis. We present the results for two approximations to the classical conformal block: up 
to the q 12 terms and up to the g 16 terms in the expansion (4.12). The results are shown on 
Fig. 2 and Fig. 3. They provide an excellent numerical verification of the classical bootstrap. 



-0.5 -0.25 

-2x10" 



S-U 



-4x10" 



tr- 



ies -Q.5' 0\75 



A- 



\ x 



4x1DJ 1 
I I 
I I 

2x10l" 
I I 
I I 
I I 



1.25 \1.5 -0.5 \-0-2,5 ~H 
\ / I ! 
Vx10( 1 l 

S-t 

-4x10"] 1 
I 



0.25 0.5 0.75 



q approximation 
-q 16 approximation 



JX, 



' \ X 



11.25 \1.5 



Fig.2: Relative deviation from the classical bootstrap equations along the real axis near x = 0, 1 singularites. 



7 Conclusions 

The numerical tests presented in this paper provide a convincing evidence that 

• the classical limit of the conformal block exists yielding a consistent definition of the 

classical conformal block; 
6 In order to simplify the formula an appropriate constant was chosen in (2.22). 
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-0.000015 



5x10' 




5x10' 



1 x10 ! 



,--0.00002 



-5x10" 



!-6 



S-U 



0.000025 A, 



S-t 



q 12 approximation 
q 16 approximation 



'SU 



-0.00001 



Fig. 3: Relative deviation from the classical bootstrap equations along the real axis near x = oo singularity. 

• the geodesic length functions of the hyperbolic geometry on the 4-punctured Riemann 
sphere can be calculated by the saddle point equation involving the classical conformal 
block and the 3-point classical Liouville action; 

• the classical Liouville action on the 4-punctured Riemann sphere can be calculated in 
terms of the geodesic length function, the classical conformal block and the 3-point 
classical Liouville action. 

The statements above were derived by heuristic field theoretical arguments within the path 
integral representation of the quantum Liouville theory and should be regarded as well mo- 
tivated conjectures. It might be surprising that they are so well supported by numerical 
calculations. Still a challenging problem is to provide rigorous mathematical proofs for them. 
Any attempt in this direction seem to require a better understanding of the classical con- 
formal block itself and in particular a direct way to calculate it. The main motivation for 
this line of research is the long standing problem of the uniformization of the 4-punctured 
Riemann sphere. Indeed if the factorization of the classical Liouville action holds and the 
classical conformal block and the geodesic length functions are available one can calculate 
the 4-point Liouville action and, by the Polyakov conjecture, the accessory parameter of the 
appropriate Fuchsian equation. 

The efficiency of numerical calculations based on the (/-expansion used in this paper are 
completely satisfactory. As an example we present on Fig. 4 the contour plots of the geodesic 
lengthes in different channels in the case of the hyperbolic metric with the constant scalar 
curvature R = — 1, i.e. the lines of constant value of the length Frenkef-Nielssen coordinate 
on the moduli space of the 4-puncture sphere as a function of the Koba-Nielsen coordinate 
x. Let us note that also the twist Frenkel-Nielssen coordinate can be calculated. Indeed it 
was proved in [24, 25] that the second derivative of the classical Liouville action with respect 
to the location of the punctures, d^d^S^ 1 ' , gives the Weil-Petersson metric on the moduli 
space. With the 4-point classical Liouville action at hand one can calculate the Weil-Petersson 
metric and the twist coordinate along each line of constant length. 
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Appendix 



The coefficients R'^ 



A 3 A 2 
A 4 Ai 



A 3 A 2 
A 4 Ai 



of the x-expansion can be written as 

A A A 



A 3 
A 4 



prs 



A 2 
Ai 



Al s = - 



dc rs (A) 



24(T rs (A) 2 -l) 



^— A rs (T rs (A)) - i _ g2 - — 2 _ l)Trs ( A) 2 A rs(Trs(A)), 



r>rs 
A 



A 2 
Ai 



(T„(A),Ai + A 2 ,Ai-A 2 ). 
For arbitrary positive integers m, n the function A rnn (a 2 ) is defined by the relations 



^mn(a 2 ) 



■I n n 

\k=l-ml=l-n/ KOt 



(A.l) 



a 



where the prime on the symbol of products means that the factors (k, I) = (0, 0) and (k, I) = 
(m,n) must be omitted. A mn (a 2 ) can be rewritten in the form 

m— In— 1 a m—l -, n— 1 



A m n(oL ) 



-. lib 1 l( _L 

inn 



n 



i 



n 



a 



fc=l «=l v ; fc=l v M=i v ; 

For arbitrary positive integers m, n the function P mri (a 2 , A, 8) is defined by the relations 

(A.2) 



a 1 + a 2 - pa + ± \ I a\ - a 2 - pa + Q 



P mn (a i ,A 1 + A 2 ,A 1 -A 2 ) = 

m— 1 n— 1 

= n n 

p=l— m g=l— n 

p+m=l mod 2 (jr+n=lmod2 

where the variables ccj are related to A, via A, = — ^ (a — + 
Pmn{a 2 , A, <5) can be expressed in the form 

4 

P mn (a,A,<5) = n^mn(«,A,S) 



fe=i 



where 



i* n (a,A,«J) = J] 



[] Q™(a,A,5)Q~™(a,A,5) 



m— l>p>0 n— l><jr>l— n 
p+m=l mod 2 g+n=lmod2 
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Pl n (a,A,S) 



11 Q°' q (a,A,S) if mis odd 

n-l>q>0 
q+n=l mod 2 

1 otherwise 



11 Q p '°(a,A,(5) if n is odd 



ra-l>p>0 
p+m=l mod 2 



otherwise 



and 
Q P 'V,A,<5) 



5 if m and n are odd 
1 otherwise 



q 2 -4 



cr 



+ (p z - 4)a z + 2(4 - pq) - 8A + S 



The coefficients 
written as 



A 3 A 2 
A 4 Ai 



Rl. 



that appear in the g-expansion of the conformal block can be 



A 3 A 2 
A 4 Ai 



j^rs prs 



A, 
A 4 



1\ _ A (u2 



A™ = A rs (-b z ) = A rs (b z ), 



A 2 
Ai 



P r ,(-6 2 ,A 1 + A 2 ,A 1 -A 2 ), 



where the functions A rs , and P rs are defined by (A.l) and (A. 2) respectively. 
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